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I. INTRODUCTION 

According to a popular viewpoint the cuprate super- 
conductors in their underdoped regime are fundamen- 
tally different from the BCS ones due to the existence of 
wide temperature range around T c dominated by phase 
fluctuations. In this range the order parameter ampli- 
tude is well defined but the global order has not yet 
emerged. The quantitative measure of this difference is 
the ratio Q = 2T c /i:p s {Q) where p s (0) is zero tempera- 
ture phase stiffness. This ratio determines how close the 
transition is to being mean-field-like. In BCS supercon- 
ductors Q <C 1, in the underdoped cuprates Q ~ 1 [|. 
Experimental evidence in favor of strong phase fluctua- 
tions comes from measurements of diamagnetic suscep- 
tibility and Nernst effect above T c [3 and from tempera- 
ture dependence of the thermal expansion coefficient Q. 
The analysis of temperature dependence of magnetiza- 
tion and London penetration depth for the high quality 
underdoped BiSCO crystals show that the superconduct- 
ing transition itself becomes closer to two-dimensional 
Berezinskii-Kosterlitz-Thouless (BKT) one Q. At last, 
there is a material (x=l/8 LSCO) where the coupling 
between the CuO planes is so weak that the transition is 
real BKT @. 

The problem of influence of the phase fluctuations on 
the spectral function is a long standing one; its impor- 
tance being substanciated by the fact that the spec- 
tral function measured by Angle Resolved Photoemis- 
sion (ARPES) serves as one of the main probes of the 
strong correlations in the cuprates. The experiments 
show that in the underdoped regime many features of 
the spectral function below T c survive above T c though 
in somewhat modified form. According to ARPES, the 
excitation spectrum above T c gradually softens and loses 
its characteristic node-centered conical shape so that the 
nodal points gradually broaden into arcs [6], 7], 9]. The 
question is whether the appearance of these arcs (or even 
of the Fermi pockets) is due to superconducting fluctua- 
tions, as suggested, for instance, in Q. To have a consis- 
tent understanding of the underdoped regime it is vital to 
obtain a correct detailed picture of the quasiparticle spec- 
tral function and establish its temperature dependence. 



The first step in this direction is to take into account su- 
perconducting fluctuations. This is possible to do even 
without full microscopic theory. The problem has been 
studied by a number of authors ( [13| , [14j , [15| ) , but, as 
we argue in this paper, the approach taken is unreliable 
being based on uncontrolled approximations. 



II. FORMULATION OF THE MODEL 

We consider a two-dimensional metal with a strong 
superconducting pairing in the state where the order pa- 
rameter amplitude is fixed, but the phase fluctuations 
are strong. Our calculations allow for the SC (super- 
conducting) gap to have nodes on the Fermi surface. For 
instance, for ci-wave SC the mean field quasiparticle spec- 
trum is given by 



E 2 = e 2 {k) + A 2 [cos(ak x ) - cos(ak y )f 



(1) 



We will assume that vja A (v is the Fermi velocity) 
and approximate the spectrum close to the node as 



E z 



l k 2 + 2A 2 sin 2 ( g a/2) 



(2) 



where k is the wave vector component perpendicular to 
the Fermi surface and q is the one parallel to it. We take 
the Fermi surface at the node for a flat one. 

In the mean field approximation fluctuations of the or- 
der parameter A are ignored, and the resulting Green 
function takes a familiar BSC form, Gfc(u;) — (u + 



e k )/(u; 2 



A (fc)). Following our original assump- 



tion we will neglect the amplitude fluctuations of the 
order parameter taking into account only phase fluc- 
tuations, A(x,i) = Ae 1 ^'*). At finite temperature T 
the long wavelength fluctuations are essentially classical 
(time independent). It is crucial for our consideration 
that the superconducting fluctuations are space isotropic 
and the one-particle Green function close to the node 
is strongly anisotropic with the parameter of anisotropy 
A/ep, where A is maximal gap and ep is the Fermi 
energy. Then when one considers a propagating quasi- 
particle, the fluctuations affect mostly its wave vector 
component perpendicular to the Fermi surface. The par- 
allel component can be considered as conserved. The 
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above considerations allow us to consider one dimensional 
fcrmions at a given Matsubara frequency to n described by 
the Lagrangian 



exponential asymptotics. The length scale a ~ (v/ep) is 
the short distance cut-off. Below the BKT transition the 
bulk Hamiltonian is Gaussian: 



Z = Xu [-iw n U-ivd a! o*+A*++A*a } X . + F . (3) H Wfc [0;T<T BKT ] = ^ J°° dy [( 47 rd) 2 H 2 + (d y ^) 2 ] ,(7) 



In the last equation the first term is a standard pairing 
Lagrangian written in terms of Nambu spinors, \lo = 
(V'tiWn ? V'4.,-ii;„) T ') arL d cr* are Pauli matrices. The second 
term gives the action for the classical phase fluctuations 
in the form 



Fa, 



2T 



dxdy [{d x <t>) 2 + (dyt) 2 



(4) 



Here the inverse temperature prefactor T^ 1 results from 
the integration over imaginary time. The discrete sym- 
metry of the lattice which includes C4, the group of 
in-plane rotations by 7r/2 radians guarantees that the 
quadratic action, Eq. ((4]) is isotropic. 

Remarkably, it is possible to rewrite the model, 
Eqs. ([3]), ((H) as an effective quantum impurity model 
of Caldeira-Leggett type [Hj]. To this end, we intro- 
duce new variables as follows, [V'ti V'j] — [*t' — *^±]) an d 

[V4>VmJ = [i*t'*|]- Then, the Lagrangian Eq. (0, © 
after the analytic continuation, iui n — > u + iO is equivalent 
to the Hamiltonian 



H eff =v-h(u) + iO)r 3 + v' 1 A(q)[r + d+W + T~e 
+ Hbuik [4>] 7 



10(0)1 

(5) 



where r a is the short hand notation for the fermionic 
bilinears:r a = ^> + a a ^ . In this setting coordinate x plays 
the role of Matsubara time. It is dual to the momen- 
tum component fcii parallel to the Fermi velocity at the 
point of observation. Since at A/ep the electron 
momentum parallel to the Fermi surface is conserved, 
the fermionic field 'I'-)-, \&t depends only on one coordi- 
nate x, though the phase field depends on both. For 
convenience we assign to coordinate y — 0. Since the 
propagators of </>-fields are not supposed to be modified 
by insersions of fermionic loops, which would lead to over- 
counting, the fermionic number must be treated as con- 
served ipeiptr = !■ Then the r-operators become compo- 
nents of spin S=l/2. The Hamiltonian H{, u ik describes 
the phase fluctuations. Their two point correlation func- 
tion is 



a 10(xi) -10(x 2 ) 



2./ 



F 



G(?)) ' 



(0) 



where d = T/8Tbkt is the order parameter scaling di- 
mension, £(T) is the correlation length and Tbkt = 
irp s /2 is the Berezinskii-Kosterlitz-Thouless transition 
temperature. The function F(y) is such that F(y <c 
1) = y ~ 2d and F(y > 1) ~ K (y) (more detailed infor- 
mation about this function can be extracted from 10] ). 
Hence below the transition where £ = 00 the function ([5)) 
decays as a power law and above the transition it has the 



where n is the momentum density conjugated to the field 
cj>, with an equal time commutator 



P(i/i),#wO] 



-i<%! - y 2 ) 



(8) 



By construction integration over the momentum in the 
quantum action of our one-dimensional model, Eqs. (J7J, 
(jHJ produces the free energy of classical two-dimensional 
thermal fluctuations, Eq. Above the transition the 
Hamiltonian Eq. ([7]) must include effects of vortices which 
generate nonlinear terms. The rigorous discussion is pos- 
sible in the case T < Tbkt! later we will present some 
extrapolation for temperatures above Tbkt- We note in 
passing that at T < Tbkt model (|5l7p is equivalent to 
the anisotropic Kondo model in the imaginary magnetic 
field h = iui — 0. We notice this equivalency though we 
have not been able to make use of it in our calculations. 

The setting of the problem as given by Eqs. (J3|), ((H) is 
not different from the one in [l3[ , [l4j , . However, in 
these previous attempts to solve the problem the authors 
used the gauge transformation method which we claim is 
inadequate. The bosonic exponent present in Eqs. ([3]). 
(O has been absorbed into the definition of the fermionic 
field. As a result the term ifj+ipo-dx^O) appeared in the 
Hamiltonian. The problem with this approach is that the 
path integral expression for the electron Green's function 
is dominated by the field configurations with large <f> gra- 
dients. The same effect appears when one attempts to 
develop a perturbation theory in d x <p: each diagram di- 
verges at small distances. This fact has been overlooked. 
Here the equivalence of the current problem to the Kondo 
model is helpful since as is well known the latter cannot 
be treated by the methods employed in [li] . fTE!] . 



III. SOLUTION BY PERTURBATION THEORY 

IN A 

In the present section we calculate the Green func- 
tion by the perturbative expansion in A for the model 
Eqs. ((3]), ((!]). It has a crucial advantage of being free of 
ultraviolet divergencies. The infrared divergencies are re- 
moved at finite external frequency and momentum (as we 
will show below, the integrals diverge only at uj+vk = 0). 
The infrared behavior is controlled by the long wave- 
length fluctuations of the order parameter which justifies 
the use of the effective low energy action for the phase 
fluctuations Eq. (g]). The above arguments indicate that 
we can study the model Eqs. ((3]), ((H in perturbation the- 
ory expanding the Green function in powers of A by ac- 
counting for most infrared singular contributions at each 
order. 
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Below for definiteness we consider positive energies, 
uj > 0. In other words we study the "particle" part 
of the spectrum. For negative frequencies the spectral 
function can be obtained by particle-hole transformation, 
uj — » — uj, and k — » —k. Most of the spectral weight is ex- 
pected to be found close to the particle mass shell, k = uj. 
In general, the Green function G£,(fc) is peaked for wave 
vectors close to the particle (hole) mass shell, k ~ ±uj. 
Here the superscript r(l) designates right (left) moving 
particles. In the former case the most singular contribu- 
tions can be resummed in a usual fashion, by introducing 
the self energy E w (fc), and writing the Green function in 
the standard form, G LJ {k) = [uj — k — £ w (fc)] _1 . The self 
energy can be written as (we set v = 1 and drop the 
explicit g-dependence of A) 



A 2 



{uj + k) 



1— 2d 



(9) 



C {d) + Y,Cn (d, 



n=l 



2uj 



A 2 



{uj + k) 



2-2d 



The coefficients C n {d, p) have singularities only at p — 
which means that the self energy is a smooth function of 
its arguments in the broad vicinity of the particle mass 
shell, k ~ uj making the resummation scheme justified. 
The direct calculation yields: 



C {d,x) 



r(l-2d), Ci{d, 1/2) «de' 



-2l7rd 



(10) 



These considerations and the fact that away from x = 
0{uj = — k) C n {d, p) vanish at d — > allows one to ap- 
proximate the Green function as 



A 2 {q)a 2d „ . ' 
UJ — k — — — — (1 — ma) 



{uj + k) 1 



-2d 1 



(11) 



It follows from equation (fTTj) that the quasi-particle dis- 
persion relation, E{k) as determined by the relation 
E(k) = k + RcE(£'(fc), k) differs from the mean field 
BSC result. In particular, the spectral gap identical to 
the Kondo scale in the magnetic impurity problem. The 
latter one is given by 



T K = A{q) [A{q)a 



d/(l-d) 



(12) 



The gap magnitude for a given q is suppressed with re- 
spect to its mean field value A(q). In addition, the spec- 
tral line acquires a finite width, proportional to the spec- 
trum renormalization T{k) w ■nd{E{k) — k) at d <C 1. 

We now turn to a discussion of the hole-like region 
k ~ — uj. The self energy is not a useful quantity in 
this case as it diverges at k = —uj and becomes more 
singular with increasing order of the perturbation the- 
ory. We therefore sum the most singular contributions 
to the Green function. For the "amputated" Green func- 
tion, Gl{k) = [G„(k) -GoAk)}{G r ,Jk)]- 2 , we obtain 
the expansion in even powers of A similar to ©. The 
coefficients of this expansion are denoted as C n . The co- 
efficients with n > 1 are non-analytic at k = —uj. To 



find the momentum dependence of the Green function at 
k ~ —uj we have calculated the whole set of constants 
C n {d, 0) and summed over the corresponding series. As 
a result of interaction with the field created by the fluctu- 
ating order parameter the right moving particle is scat- 
tered off as a left moving hole, see Fig. [JJ When the 



(a) 
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FIG. 1: (Color online) Propagation of the right-moving par- 
ticle, (a) k ~ +Lu and (b) k ~ —uj. Solid straight arrowed 
(black) lines labeled with letters r (V) designates propagators 
of the right-moving particle (left-moving hole). Pairs of out- 
going (red) and incoming (blue) dashed wavy lines designates 

the insertion of Ae 1 ^ and Ae -1 * 1 phase factors with close ar- 
guments. In the panel (b) a pair of distant solid wavy lines 
depicts the two uncompensated phase factors. 



order parameter is fixed these processes give rise to a 
usual BCS-like quasi-particle spectrum. In our case the 
order parameter has different phases at different collision 
events. Our goal is to study the effect of these fluctua- 
tions to the leading order in the parameter d. For uj ~ k 
the right-moving particle propagates along much longer 
distance in between the consecutive scatterings than the 
left moving hole. This is shown in Fig. Ufa). As a re- 
sult the space arguments in the phase factors attached 
to the propagator of the hole merge. For that reason 
the integration over the short distances of propagation of 
the hole leads to effective fusion of its propagator. This 
explains why the Green function is determined by the 
self energy to the second order in A, i.e. the result of 
Eq. (fTTj) . Indeed, corrections to Eq. (fTTj) results from 
the interaction between distant dipoles and are small in 
parameter d ~ T. 

The situation at uj ~ —k is different, see Fig. [2(b). 
In this case the left-moving hole propagates much longer 
distances than the right-moving particle. As a result the 
arguments of the two outermost phase factors are sepa- 
rated by large distance which substantially modifies the 
infrared behavior of the propagator. The two uncom- 
pensated phase factors lead to a power law dependence 
in Eq. (fTT)) absent in Eq. (jlip . We now derive the result 
of Eq. (fTT|) analytically. To illustrate the calculation we 
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consider the order A 6 , (see Fig. QJb)): 

SGl(R) = (-i) 6 A 6 G l 0i _ u (yi - R)G^( yi - Xl ) 

xGo,-«(ifc - xi)G r 0tM {y 2 - x 2 )G\_ u) {-x- 2 ) 
x ( e -U(R) e U(yi) e -U(xi) e my2) e i<P(x2) e -mo)^ _ q 3 ) 



*0,±u ( x ) 



Here the bare retarded Green functions ±iG 

6(x)e lulx . The infrared singularity in the Fourier trans- 
formation of Eq. (fT3|) at k ~ — cj accumulates at large 
distances, i? cx — \ut + k^ 1 . In terms of rescaled vari- 
ables, 211,2 = — R£i,2> J/1,2 — %i,2 — —Rf]\,2 it means that 
the leading contribution comes from 771,2 <C 1. Integrat- 
ing over negative R's we obtain 



d77id?72 



(p + m + mY~ M 



(14) 



As the integrals over 77 variables are convergent allowing 
us to replace the upper integration limit by infinity. The 
remaining integration in Eq. (|14l) yields 



tf( 6) G£(fc) 



A 6 r 2 (l - 2d)T(3 - 2d) 
2!(2w) 5 - 6 V" 2d 



(15) 



The arguments leading to Eq. (|15p can be generalized to 
obtain the most singular contribution to arbitrary order 
A 2 ™ as follows 



6 i2n) GZ(k) 



A 2n r n - 1 {l-2d)T{n-2d) 



(n - i)!(2w)2«(i-d)-yi-2d ' 
Summing all leading singularities Eq. (|16[) we obtain 



(16) 



A 2 (g)r(l -2d)e 



-lird 



ZUJ |U+ fc ( 2w )i-2d 



l-2<2 



(17) 



At d — > (|17[) reproduces the BCS Green function. 
(fTTj) . (fTT)) are the main results of the paper. These results 
remain qualitatively correct above the BKT transition, 
provided that the inverse correlation length £~ 1 (T) ex- 
ceeds the Kondo scale Tk given in Eq. (|T2|) . Since the 
BKT correlation length is exponentially large £ _1 (T) ~ 
Ao exp [— C(T/Tbkt — l) -1 / 2 ] , there is a range of tem- 
peratures and q where this condition is fulfilled. In the 
opposite limit 3> A(g) one can use Eq. (|BJ)) to cal- 
culate the leading order contribution to the self energy, 
The following formula provides an interpolation between 
T > Tbkt and T < Tbkt region: 

k,«j) = A 2 (qm/a)- 2d (l + (tfw + fc)) 2 )- 1 / 2 

(18) 

x [exp {— md + 2dsinh _1 [^(tj + fc)]} 



(l + (£(u; + fc)) 2 ) 



2\-l/21 




FIG. 2: (Color online) The spectral function at k — as 
a function of dimensionless parameters q = A(qa^/v) and 
u] = (lj£/v) as found form Eq. ([18} for a) d= 0.125 and b) d= 
0.5. 




FIG. 3: (Color on line) The dependence of the frequency w ma x 
maximizing the spectral function (|18l) on the dimensionless 
momentum q for d — 0.5 solid line (black), d = 0.25 dashed 
line (blue), and d = 0.125 dashed dotted line (red). 



In Fig. [5] we present graphically the spectral function 
Ad(fc) = — (l/7r)ImG w (fc). The quasiparticle dispersion 
can be identified as the energy w max where the spectral 
function is at its maximum. The dispersion relation ob- 
tained in this way is depicted in Fig. [3] for different val- 
ues of parameter d. Below T < Tbkt where = 
one should use Eq. (|T7)) instead of Eq. ([TBI close to the 
singularity line u) — — fc. 



IV. CONCLUSIONS 

In summary, we have studied the electron Green func- 
tion (and the related spectral function) in the regime 
of strong superconducting fluctuations. As is evident 
from Figs. [5] and [31 these fluctuations affect the nor- 
mal state dispersion in such a way that the maximum 
of the spectral function is shifted down in energy in com- 
parison with its mean field value. Naturally, the effect 
is more pronounced for larger temperatures. This in- 
deed may create the impression that the system develops 
a "Fermi arc". We hope that the advances in the ex- 
perimental techniques will allow for a detailed compar- 
ison with the present theory. We argue that the quasi- 
classical approximation employed in the previous publi- 
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cations [13 > [3 > [3 cannot provide a quantitative infor- 
mation about the spectral density along the arcs. This 
approximation is justified only if Green function changes 
on the scale smaller than the variation scale of the pairing 
potential. In general the former is set by the particle's 
mass. As the quasi-particles become massless at the node 
the quasi-classics is not justified when A(k)/v < £ _1 (T). 
More specificall y th e inverse square root singularity re- 
ported in Ref. IS] is an artifact of the quasi-classical 
approximation and is in fact smeared. The typical scale 
of the spectral function is temperature dependent and 
can be estimated as ~ <i(£/v) A 2 (?;/£) close to the node. 

Our approach is easely generalizable for other problems 
where quasiparticles coexist with critical or almost crit- 



ical collective excitations such as magnetic fluctuations 
at the onset of antiferromagnetic order. 
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